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Abstract 

We consider the effect of different unitary noise meclianisms on tlie evolution of a 
quantum walk (QW) on a linear chain with a generic coin operation: (i) bit-flip 
channel noise, restricted to the coin subspace of the QW, and (ii) topological noise 
caused by randomly broken links in the linear chain. Similarities and differences in 
the respective decoherent dynamics of the walker as a function of the probability 
per unit time of a decoherent event taking place are discussed. 
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1 Introduction 



Quantum walks (QW) [T] are a generalization of random walks to the quan- 
tum mechanical regime. Due to quantum interference effects, the QW spreads 
faster than its classical counterpart [2||3] . This has motivated several optimal 
quantum search algorithms based on discrete-time |^||5] and continuous-time 
versions ^ of the QW. Other recent algorithmic applications of quantum walks 
include the best algorithm for the problem of element distinctness {i.e. deter- 
mining if all elements in a set are distinct or not) [7J and an algorithm for fast 
evaluation of NAND trees [8]. Furthermore, a continuous time quantum walk 
has been shown to traverse specific binary tree-shaped networks exponentially 
faster than a classical random walk [QllTO] . Even though the significance and 
generality of some of these results is still under debate it is clear that 
QW's are a highly useful concept for quantum computation. 
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Among the many suggestions made for the physical reahzation of the QW, 
some are now implemented at the proof-of-principle stage [T2p3] . As in any 
implementation of a quantum device, the problem of decoherence due to cou- 
pling to the environment or to imperfect gate operations, is a mayor obstacle 
to be dealt with before useful computations can be accomplished using QWs. 
Quantum error correction protocols require a considerable overhead in quan- 
tum resources, so it is important to characterize the effects of decoherence on 
the dynamics and perhaps even use controlled decoherence to achieve specific 
purposes In recent years, several studies on decoherent quantum walks 

in one [TB|17|18II19|2U] or more dimensions pT|22y23] have appeared. Ref. [21] 
provides a recent review on the subject. 

The effects of decoherence on the quantum dynamics are frequently studied 
by applying some quantum operation on the system with probability p per 
unit time. When the quantum operation is restricted to the coin subspace, 
considerable progress can be made using analytical techniques. In other cases, 
path counting-methods have proved to be useful tools, but frequently one must 
resort to numerical simulation to explore the resulting dynamics in detail. It is 
a common feature of decoherent quantum walks that for times long compared 
to a characteristic time 1/p, the variance increases linearly, but at a higher 
rate than that implied by classical diffusion. This spreading rate is a useful 
signature of the effect of decoherence on the quantum walk. 

Usually some fixed unitary coin operation is used to drive the quantum walk. 
However, the impact of decoherence may be different for different coin op- 
erations. In this work, a decoherent quantum-walk on the line driven by a 
generalized coin operation is analyzed. Two different noise mechanisms are 
considered: (i) a coin-flip operation applied with probability p per unit time 
{i.e. bit-flip channel noise) and (ii) at a given time a link on the line is open 
with probability p. This model (broken-link model) was introduced in [20] as a 
way to mimic the effects of thermal noise in some experimental situations and 
it was later generalized to higher dimensions . This work is organized as 

follows: in section [2] the formalism of the coherent quantum walk is introduced 
and, in this context, the method of quantum operations is reviewed. In sec- 
tion [3] the two specific noise mechanisms mentioned above are considered and 
their effect on the spreading rate of the quantum walk is discussed. Finally, in 
section H] we present our conclusions. 



2 Generalized discrete-time quantum walk on the line 

The discrete-time quantum walk on the line is a quantum analog of the classi- 
cal random walk where the random choice is replaced by a unitary operation 
in the abstract "coin" subspace, Ttc- This is a single-qubit space, spanned by 
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two ortho normal vectors, usually denoted |L)}. In the one dimensional 
case, the motion takes place in a position subspace. Tip, spanned by an or- 
thonormal set of position eigenstates {\x)} with x an integer associated to 
discrete positions on a hne. The Hilbert space is then Ti, = Tip ® Tic and a 
generic state for the walker is of the form 

oo 

l^)= E \x)0ia,\R) + K\L)), (1) 

2:=— oo 

where the amplitudes = {x, R and = {x, L satisfy the normaliza- 
tion condition J2x + l^xP = 1- 

A step of the walk is described by the unitary operation 

U = So-{Ip® Uc) (2) 
where Ip indicates the identity in Tip and 5*0 is defined by 

So = Y. (1^ + ® + 1^ - ® \L){L\) ■ (3) 

X 

This unitary operator conditionally shifts the position by one step. The coin 
operation, Uc, can be any suitable unitary operation in Tic- 

2.1 Generalized coin operation 



A parametrization which represents a general unitary operation in Tic requires 
three real parameters. For our purposes, a single-parameter coin operation that 
allows us to explore the response of the system to different environments will 
suffice. We consider coin operations 



c 



cos 9 sin 9 
sin 9 — cos 9 



(4) 



parametrized in terms of the real angular parameter 9 G [— 7r/2, 7r/2]. For 
9 = 7i/4, the QW reduces to the standard Hadamard walk on a line. 

Open quantum systems are best described in terms of density operators. In the 
absence of noise, an initial state po evolves, after t iterations, to pt = U^po{Wy. 
The probability distribution for finding the walker at site x at time t is 
P{x,t) = tr{pt\x){x\), where tr(-) represents a trace operation. It is a well 
established fact [2126 ] that the variance of this distribution increases quadrat- 
ically with time, while in the classical case the increase is only linear. The 
quadratic increase is directly related to quantum interference effects and is 
eventually lost in the presence of decoherence 
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Since the QW has a constant step size, its dynamics is best described in the 
Fourier-transformed space [2^. This subspace, Tiki is spanned by the Fourier 
transformed kets \k) = J2x^^^^\^)j where the real wavenumber k is restricted 
to [— 7r,7r]. In this representation the density operator is 

P= 7^ -7^\k){k'\®Xkk' (5) 



TT 2tT J-tt In 

where Xfefc' = \Xklkx.k'\- In this work we consider pure, locahzed initial states, 
and, without loss of generality, start the walker at the origin, x = 0. The ini- 
tial state is, therefore, of the form po = |^o)(^o| with |\l/o) = |^|^) ® Ixo), 
where |xo) is an arbitrary initial coin state, with components (Sfc, h^)^ = (ao, ^o)"^ 
satisfying |aoP + |&oP = 1- We use (a, b)^ to indicate_a two-component column 
vector and the amplitudes in Fourier space are {ak,bk)^ = J2x^~^^^{C'x,bx)^ ■ 

The shift operator is diagonal in Hk since So\k,R) = e''^^\k,R) and 

So\k, L) = e^^\k, L). Therefore the evolution operator, eq. ([2]), is diagonal in k 

and acts non-trivially in the coin subspace, U{\k) ® \xk)) = \k) ®Uk\Xk), with 



e cos 9 e ^'^ sin 9 



(6) 



2.2 Quantum operations 



Long-time effects of different kinds of noise on the position distribution of 
the walker can be obtained using quantum operations, as described in detail 
in Ref. [17]. Here, we review the essential aspects of this method and intro- 
duce the appropriate notation. A trace-preserving quantum operation [28] is 
described by a set of Kraus operators {^n}, n = 1, 2, . . . assumed to satisfy 

N 

E ^i^n = I. (7) 
n=l 

Considerable progress in an analytical description of the dynamics is possible 
when the operators An are restricted to the coin subspace He- In this case, 
Xkk' — ^ J2n ^nXkk' and, after a single step in the quantum walk, the coin 
state becomes Xkk' J2nUkAnXkk'Al^ul,. After t such steps are taken, the 
density operator is of the form ([5]), with 

Xkk'it) = E UkAn, . . . UkA^^xoAlul . . . Alul,. 

ni,n2...nt 
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This transformation is best described in terms of a superoperator Ckk' defined 
by 

N 

Ckk' [Xkk'] = UkAnXkk'Alul, (8) 

n=l 

SO that Xkk'(t) = i^\k' [Xo]- Notice that the diagonal instance of this superop- 
erator, Ck = Ckk, is trace preserving. 

We are interested in the effects of decoherence on the long-time position proba- 
bility distribution, P{x, t), of the QW. In what follows, unless otherwise stated, 
the trace operation acts only in the coin degree of freedom. The moments of 
the position distribution can be obtained from the above expressions [TT], 



X 

=E^"r^r ?^^''^'-''^tr[xkk'{t)] 

^ dk'6^"^\k-k')tr[xkk'it)], 

where the sum over positions has been evaluated in terms of derivatives S^"^^ 
of the delta function. After integration by parts the first two moments, which 
determine the variance cr^, are given by 



(x) 



and 



X 



j'=i 

+ x:^-{z^r'((4[xo])z)} 



where the fact that Ck is trace-preserving has been used. 



(9) 



(10) 



In order to find a clean expression for C^ [x], the reduced density operator 
X is parametrized as a linear combination of Pauli matrices and the identity, 
X = I]j=o'^iCri5 where ctq = I is the 2x2 identity and cxi = X, c72 = Y and 
(T3 = Z are the usual Pauli matrices. Since tr{x) = 2ro = 1, ro = 1/2 is 
fixed and Ck does not affect this component, the state of the walker can be 
parametrized by the column vector R = (ri, r2, rs)"^. Then, the action of the 
superoperator Ck is described by a 3 x 3 matrix Mk, 



R' = MkR. 



:iii 
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For a specific set of Kraus operators {^n}, tlie matrix representation Mk can 
be obtained from eq. ([8]). We perform an explicit calculation of this kind in 
the following section. 



Moments of the position distribution 

In the absence of decoherent events, the (long-time) first and second moments 
are proportional to t and t^, respectively. In the presence of decoherent noise, 
the moments of the position distribution may be expressed in terms of Mk as 
follows. 

The first moment, eq. ([2]), requires the evaluation of tr |Z£;i[x]| in the Pauli 
representation. Since tr |Z£|.[x]| = 2r'^, the first moment reduces to 

-Y.Mi{r,,r,,r,r. (12) 

Provided the eigenvalues Aj of satisfy |Aj| < 1, at long times the sum may 
be approximated by a geometric series which sum is given by the constant, 

— Gfc(ri,r2,r3)^. (13) 

We have defined the operator 

Gu = {I-Mk)-^Mu (14) 

assuming that I — M^, is invertible (/ is the 3x3 identity matrix). Note that, 
due to the scalar product with the row vector (0,0,2), only the third row of 
Gk is relevant. We emphasize that this expression for the first moment of the 
decoherent evolution is valid only in the limit of very large times. 

A similar expression for the second moment can also be obtained, although the 
details are more involved. The diagonal term {j' = j) in eq. ffTOl) is readily eval- 
uated, and represents the "classical-like" linear contribution to the variance. 
The remaining terms represent quantum effects and and may be regrouped in 
the form 



i^') = ^ + 1^ E E {Z^r" [Z (4 [xo]) + (4 [xo]) Z] } . (15) 



dk * 
j=i j'=i 

This expression may be written in the 4x4 Pauli representation as 

t j-i 

EE 

I i=i j'=i 



Xo, (16) 
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where v = (0, 0, 0, 2) and the matrices Zl and represent the action of 
the Pauh operator Z on the left or right, respectively, i.e. Z^O = ZO and 
OZr = OZ. The explicit form of these matrices is 



A 
-i 

i 

1 Oy 



Zr = 



A 
i 

-i 

1 Oy 



(17) 



so Zr + Zl is sparse and real. The second moment should not depend on 
the initial conditions. As shown in Ref. [17], this can be made apparent by 
separating the first component of the initial coin state 



Xo = (l/2,0,0,0)^+(0,ri,r2,r3)^ 

and noting that 

V ■ Lr\ZL + ^ii)4'(0, ri, r2, r^f = v ■ (1/2, 0, 0, 0)^ 



:i8) 



0. 



After taking into account the action of Zl + Zr on the first term of eq. f[T8|) . 
the expression for the second moment may be further simplified to 



{x')=t + v. r|^EE^r'(o, 0,0,1)^ 

=t + (0,0,2). /^^Ex:Mr'(o,o,i)^ 

where the first trivial component has now been omitted. The double series 
may be summed as before, provided the eigenvalues of Mk satisfy the restric- 
tion I Ail < 1. Then, in the limit of very large times, the second moment is 
approximated by 



(x') (0,0,2) 



^ dk 



tl-il-Mk)-^ Gk (0,0,1) 



(19) 



Thus, the variance of the QW is determined by the element in the third row 
and column of Gk- 



Spreading rate 

As we have seen, in the presence of decoherence the first moment is a constant 
and the second one increases linearly with time. When the rate of decoherent 
events is small this increase can be considerably faster than the corresponding 
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classical diffusion rate [20] . We define the spreading rate of a quantum walker 

as 

D, = \im-^. (20) 
^ t»i dt ^ ' 

If the variance o"^ is associated to an ensemble of walkers, the corresponding 

classical quantity is the diffusion coefficient of a random walk, D = 1. Since in 

this work t is a discrete variable, the use of discrete derivatives in the above 

definition is implied. As mentioned before, the first moment is asymptotically 

constant and does not contribute to Dq. A simple explicit expression for the 

spreading coefficient can be obtained from eq. f|T9l) . 

= 1 + 26-3,3. (21) 

The upper bar indicates that an average in fc— space has been performed, i.e. 
G = J^T, ^Gk- Thus, the term G3 3 is responsible for the faster spreading rate 
of a decoherent quantum walk relative to a classical walk. 

For any particular kind of noise, eq. ([8]) allows one to evaluate eq. (12T|) and 
obtain Dq in terms of the noise rate and coin operation. In Ref. [T7| this was 
done for the particular case of a Hadamard walk {i.e. 9 = 7r/4 in eq. (^) 
and for a set of Kraus operators representing a partial measurement of the 
coin state with probability p per time-step. In the following, we apply this 
formalism to other sources of decoherent events and arbitrary coin operations 
of the form (HI). 



3 Decoherent quantum walk 



3.1 Bit-flip channel 



Consider the particular set of Kraus operators which flip the coin state with 
probability p per time step {p G [0, 1]), 

Ao^v^X, Ai = yr^/. (22) 

This particular kind of noise is usually known as bit- flip channel noise [28j. 
Note that these operators satisfy eq. ([7]) as required from a trace-preserving 
quantum operation. For the bit-flip channel, eq. ([8]) reduces to 

X' = p f/fcX X Xf/i + (1 - p) Uk X Ul (23) 

with Uk given by eq. ([6]). There are two parameters in this model; the probabil- 
ity per unit time that the coin state is inverted, p, and the angle 6 which deter- 
mines the coin operation. The purity of the state p, defined as Il{t) = tr{pf), 
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Fig. 1. (color online) Purity vs. time for a Hadamard walk {0 = tt/4:) with bit-flip 
noise. Several noise rates are shown: p = 0.01 (full line), p = 0.03 (dashed), p = 0.05 
(dash-dot), p = 0.10 (dash-double dot) and p = 0.50 (dot-double dash). The asymp- 
totic power-law dependence ~ is shown as a dotted line to guide the eye. 

can be used as simple indicator of the impact that a given noise rate has on the 
dynamics. Fig. [T] shows the purity as a function of time for several noise rates 
p < 0.50. For values of p > 0.50 the purity is identical to that of 1 — p. For 
all noise levels and times long compared to 1/p, the purity decays according 
to a power law, namely 11 ~ This decay rate was to be expected since, 

for t 3> the density operator reduced to the coin subspace, pc = trc{p), 
has evolved to a minimum information state, i.e. pc //2. In this regime, 
the position distribution is gaussian-like, peaked at x = 0, with a characteris- 
tic spread cT{p,t) = JDgt, as implied by definition (!20l) . A simple calculation 



shows that this leads to a purity decay 11 ~ 1/ \J Dqt. 

We now proceed to obtain an explicit expression for the spreading coefficient 
Dq{p, 9). After expressing x in the Pauli representation introduced in the last 
section and using the properties of the Pauli matrices, we obtain from eq. flTTl) 
and eq. (1251) the expression for M^, 



v 



cos26'cos2fc gsin2A; gsin2^cos2A; 
■ cos2^^sin2A; —qcos2k qsm26sin2k 
sin 2^ gcos2^ 



(24) 



where q = 1 — 2p. The operator / — Mj. may be inverted provided 
det(J-Mfc) = (1 -g^) [1 + cos2^cos2A;] =^ 0. 



(25) 



This condition is satisfied provided p > 0, 9 ^ and 9 ^ ±7r/2. We shall 
discuss these special coin operations later. 
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Fig. 2. Dependence of the spreading rate Dq on the probabihty p, from eq. (j26p . 
Three coin operations are shown: = 7r/4 (Hadamard coin, full line), 9 = 7r/30 
(dashed line) and = 0, eq. (f28l) . (dash-dot line). 



Our main interest is to find the dependence of the spreading coefficient on 
the parameters p and 9. According to eq. fl2Tl) . the relevant information is 
contained in Gs^s which equals 



G 



3,3 



q dk {1 + q cos 29) cos 2k + q + cos 29 



1 - g2 2 



1 - 



in 



1 + cos26'cos2/c 



1 - I sin 2^1 
cos 20 



Thus, the spreading rate for the case of bit-flip channel noise is. 



D,{9,p) = -^^ + 
I — q^ 



2q 



sin 201 



cos 29 



(26) 



This expression is one of the main results of this work. One of its interesting 
features is that Dg = 1 for p = 1/2, regardless of the coin operation. This can 
be understood from the fact that p = 1/2 represents minimum information in 
eq. (!23ll and corresponds to the highest decoherence rate. It also explains the 
fact, apparent in Fig. [H that lower p values correspond to lower purity at long 
times due to their faster spreading rates. For p = and for p = 1, eq. ( I26l) does 
not apply and in fact diverges. This is to be expected on physical grounds, since 
in these cases the quantum dynamics is coherent and the variance increases 
quadratically with time. The case p = 1, which corresponds to a coherent 
quantum walk with a modified coin operation Uq = UcX, has been considered 
in detail in 



For a Hadamard walk, 9 = 7r/4, the second term in the expression for the 
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9 

Fig. 3. Spreading rate, for coin operations of the form Q for noise from different 
sources. Left panel: Bit-flip channel noise for noise levels p = 0.01,0.02 and 0.10. 
The full lines correspond to the asymptotic {t — > oo) expression for Dg, eq. (p6]) . 
The dots are estimates of Dq obtained from eq. (I20p through linear regression, after 
following the time evolution of an ensemble of 1000 walkers for 1000 time steps. The 
discrepancy as approaches ib7r/2 is a consequence of this finite time, as commented 
in the text. Right panel: Topological noise (broken links) for the same noise levels 
as in the right panel. In this case Dq was obtained numerically iterating the map 
described in Ref. [20J for 1000 steps for an ensemble of 1000 walkers. In both cases, 
the maximum spreading rates are consistent with a (1 —p)/p dependence. 

spreading rate is null, so (126|1 simplifies to, 

which for low bit-flip rates, p <^ 1, is of the order of l/(2p). The case p = 1/2 
corresponds to a minimum spreading rate as shown in Fig. [2] and higher or 
lower values of p correspond to lower decoherence and faster dispersion rates. 

For ^ = the coin operation is diagonal and the quantum walk is equivalent 
to a classical random walk in which the walker reverses its direction of motion 
with probability p per time step. This walk has a spreading rate 

Dq = = for ^ = 0. (28) 

1 — g p 

For 6 = ±7r/2, the coin operation reduces to X and the quantum walk is 
equivalent to a classical random walk in which the direction of the walker is 
preserved with probability p. In this case, the spreading rate is 



Dq = - — ^ = for e = ±.1x12. (29) 



For small values of p, this leads to a very slow spreading rate, a factor of p^ 
smaller than the maximum rate, eq. (!28|) . 

We have checked the Q dependence in eq. (!26|) by performing an independent 
numerical simulation of the quantum walk on the line in which the coin is 
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inverted with probability p. The wavevector at time t results from 

|^(t)) = f/t...t/i|^(0)) (30) 

where each unitary operator Uj is either U with probability 1 — p, 
or U' = {Ip ® X) ■ U with probability p. For times long compared with 1/p, the 
slope of 0"^, and thus an estimate for Dq, are obtained by linear regression. An 
ensemble of walkers with the same initial condition was used in order to reduce 
fluctuations. The resulting dependence of Dg on 9 is shown for several p values 
in the left panel of Fig. [3] with good agreement with the analytical result. The 
dispersion rate tends to be underestimated by the simulation near the singular 
points 6 = ±7r/2 due to the fact that the numerical estimate is based on 
finite-time evolutions. As 6 approaches these singular values, larger times are 
required to adequately describe the dynamics since the analytical expression 
corresponds to the limit t —>■ oo. We have checked that if larger times are 
considered, the difference with the theoretical rate is reduced accordingly. 

It is illustrative to compare the spreading rate obtained for bit-flip channel 
noise with existing results for other noise channels. Noise events can be distin- 
guished according to whether they are restricted or not to the coin subspace 
of the QW. Among the first kind, Brun et al. [17] obtain analytical results for 
the dispersion rate of a Hadamard quantum walk (HQW) due to partial coin 
measurements (or, equivalently, pure dephasing noise). Specifically, they con- 
sider the quantum operation defined by Aq = ^/2p\R){R\, Ai = ^/2p\L){L\ 
and A2 = y/l — 2p I. The spreading rate implied by the results of that work is 
identical with our eq. (!27|) . However, it is important to emphasize that the bit- 
flip channel noise considered here is an example of unitary noise, in the sense 
of eq. (l30l) . where no measurements are made during the evolution. Shapira et 
al. have previously considered the effect of another kind of unitary noise (also 
restricted to the coin subspace) on the HQW [18] . In that work, in addition to 
the Hadamard coin operation, a stochastic unitary operator e*"^ was applied 
at each step. The hermitic operator A = aiX + + a^Z was defined at 
each step by choosing at random the coefficients <yj{t) with variance = 2p. 
For small p, the numerical results in [18] imply that Dg l/{2p), which is 
consistent with the weak noise rate limit of eq. (127|) . 

A different situation arises when the decoherent events are not restricted to 
the coin subspace. Kendon and Tregenna have considered both numerically 
and analytically [MfTB] the case of complete measurements (both coin and 
position) performed with certain probability. They investigate two limiting 
cases. In the first, the walker takes a number of steps <^ 1/p and those results 
do not apply to the long time limit considered in this work. The other corre- 
sponds to the classical limit, g ^ 0. In this later case, their results imply a 
spreading coefficient Dg ^ 1 + q'^, while from eq. (p7|) we find Z)g ~ 1 + 2g^, 
for the case of unitary noise. 
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3.2 Topological noise from broken links 

A different Icind of noise, affecting both the coin and position subspaces, ap- 
pears when the hnks between neighboring sites in the hnear chain are broken 
at random. The broken-hnks noise model was introduced in [20] for the quan- 
tum walk on the line and generalized to more dimensions in [2T]. It assumes 
that, at a given time t, a site x in the line has one of its neighboring links open 
with probability p. A broken link remains so for a unit time step, so p is in fact 
a probability per unit time. If a link is broken, the corresponding probability 
flux is not transferred across it and is diverted to the other coin amplitude 
at the same site. This is done using, in addition to 5*0 given by eq. Q, the 
additional shift operators 

X 

= E ® + 1^ - ® ' (31) 

X 

'53=(ei^)(^|)®(|i)(o| + |o)(i|) = ^®x. 

During the evolution, each of the operators Sq.Si, S2, S3 is applied with prob- 
abilities (1 — py, p{l — p), p{l — p) and respectively, so a set of Kraus 
operators for this kind of noise is 

A, = {l-p)S, Ai = v'p(l-p) 5i, A2 = ^p(l-p)^2, A3=PS3. (32) 

This set of operators satisfies eq. ([7]), as required for a trace-preserving quan- 
tum operation. The shift operation follows the unitary coin operation Uc{0), 
defined in eq. @, so that a step in the evolution can be represented by 
p = EpE'^ with 

E={j2A^)-{Ip®Uc). 

n.=0 

Since the Kraus operators defined by eqs. fl32l) are not restricted to the coin 
subspace, the theoretical formalism described in the previous section cannot 
be applied. 

An alternative approach is to describe the evolution as the result of applying a 
sequence of uncorrelated unitary operators on the initial state, as in eq. fl5U]) . 
Now each Uj describes the unitary operation that takes place when some 
specific links in the line are broken at timestep t = j. This state-function 
approach leads to a map for the wavevector amplitudes, as shown in Ref. |20j . 
In that work it was established that for a HQW with randomly broken links 
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the spreading rate depends on p as 

D, oc (33) 
p 

with a proportionahty constant of order 1. A similar dependence has been 
obtained for higher dimensional quantum walks when the links are broken 
isotropically [2T]. Here, we are interested in how Dg depends on the choice 
of coin operation as determined by the parameter 6. This can be obtained 
numerically by following the dynamics of an ensemble of walkers to times 
t ^ 1/p for several values of 9 and calculating Dg from the cr^ vs t data, using 
eq. (!20|) . Our results, for different noise levels, are shown in the right panel 
of Fig. [31 Comparing both panels in this figure, we notice that while both 
noise models have a similar dependence on the noise rate around ^ = 0, the 
details of their behavior in this region are quite different. While the broken- 
links noise is flat, the bit-flip noise presents a discontinuity in its derivative, 
specially noticeable at weak noise rates p <^ 1. 

The special cases 6 = ±7r/2, for which Dg is very small, have been discussed in 
the previous section. For the case ^ = 0, where the coin operation reduces to 
Z, the spreading rate is fastest as was in the case of bit-flip channel noise. In 
this case the motion of the individual walkers is very sensitive to decoherent 
events and leads to large fluctuations in the values of Dg obtained through 
numerical simulation. The Hadamard walk corresponds to an intermediate 
spreading rate. For instance, the HQW spreading rates may be doubled if coin 
operations with small 9 values are used. 



4 Conclusions 

The decoherent quantum walk (QW) on the line has been investigated, both 
analytically and numerically. At long times compared to the decoherence rate, 
the QW position distribution spreads out with a variance proportional to 
the discrete time t. The proportionality constant is the spreading rate, Dg, 
which depends both on the noise rate p and on the particular noise model 
considered. The noise models can involve complete or partial measurements 
or the application of different unitary operations at each time step (unitary 
noise). In this work, the spreading rate of a generalized QW resulting from 
two kinds of unitary noise has been investigated: (i) bit-flip channel noise, 
restricted to the coin subspace of the QW and (ii) broken-links noise, which 
affects both subspaces. 

A family of coin operations generated by a single real parameter 6 has been 
considered. This family includes special cases, such as the Pauli operations 
X and Z, and the Hadamard coin, a preferred choice for the one-dimensional 
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QW. The dependence of the spreading rate on this family of coin operations 
has been investigated for both kinds of noise. In the Hadamard case, 6 = tt/A, 
our results are compared with existing results for other noise models. The 
spreading rate is highest for ^ ?s 0, which corresponds to essentially diagonal 
coin operations, close to Z. In this case, an initial coin superposition state is 
preserved throughout the evolution and, in this sense, the walker experiences 
a simple translation. For coin operations with 6 ~ ±7r/2, or close to X, in 
the coherent case {p = 0) the walker remains "locked" in the neighborhood 
of the starting point and no spreading occurs. For p > 0, a decoherent event 
"unlocks" the walker, but the spreading rate is a minimum, independently of 
the noise model. Due to these general features, the overall dependence of the 
spreading rate on the coin operation is independent of the noise model. The 
details are, however, different when noise affects the position subspace and 
when it does not. 

If the decoherent events are restricted to the coin subspace, considerable 
progress can be made using quantum operational techniques. For the QW 
with bit-flip channel noise an analytical expression for the dependence of the 
spreading rate on the coin operation and on the noise rate has been obtained. 
This expression, valid for long times, was found consistent with the results of 
finite-time numerical simulations. The spreading rate increases as 1/p for weak 
noise and, as expected, approaches the classical value Dg ^ 1 for high noise 
rates. Within this noise model, the purity 11 = tr{p^) of the density operator 
p of the QW was found to decay according to a simple power law, namely 
n ~ t"^/^. For the particular case of a Hadamard walk, the dependence of the 
spreading rate on the probability p is consistent with the results reported in 
pT] which uses partial measurements of the coin state as the source of noise. 
This suggests that the dependence on the frequency of decoherent events may 
have a generic character or at least be weakly dependent on the details of the 
quantum operation, when it is restricted to the coin subspace. 

When the noise model affects the position space, the dynamics is different 
from the previous case. For instance, important differences have been reported 
in the position distribution using frequent position measurements as a noise 
model [H] . The effects on the spreading coefficient of a different noise model, 
the broken-link model, which affects both the coin and position subspaces of 
the QW were investigated and compared to that of the bit-flip noise. The 
comparison is of interest, because both noise models correspond to unitary 
noise and involve no measurements, but one of them (bit-flip) is restricted to 
the coin subspace. A different dependence on the coin operation parameter 
9 was found, but the dependence on the noise level parameter p is similar in 
both cases, particularly at low noise levels. 

The previous considerations may be summarized by noting that the noise rate 
p is of primary importance in determining the dynamics of the decoherent QW. 



15 



In all cases considered, the spreading rate is higher than the corresponding 
classical rate by a factor of the order oil/ p. In a previous work [20] we argued 
that this is an indication of persistent quantum correlations. In other words, 
even for arbitrarily long times, partial quantum coherence continues to play a 
role in the dynamics of the evolution. A similar conclusion is obtained for the 
coin-flip noise introduced in the present work. 
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R.D. and G.A. acknowledge support from the Milenium Institute for Quantum 
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